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Abstract. Observations from the two STEREO-spacecraft 
give us for the first time the possibility to use stereoscopic 
methods to reconstruct the 3D solar corona. Classical stere- 
oscopy works best for solid objects with clear edges. Con- 
sequently an application of classical stereoscopic methods 
to the faint structures visible in the optically thin coronal 
plasma is by no means straight forward and several problems 
have to be treated adequately: l.)First there is the problem 
of identifying one dimensional structures -e.g. active region 
coronal loops or polar plumes- from the two individual EUV- 
images observed with STEREO/EUVI. 2.) As a next step one 
has the association problem to find corresponding structures 
in both images. This becomes more difficult as the angle 
between STEREO-A and B increases. 3.) Within the recon- 
struction problem stereoscopic methods are used to compute 
the 3D-geometry of the identified structures. Without any 
prior assumptions, e.g., regarding the footpoints of coronal 
loops, the reconstruction problem has not one unique solu- 
tion. 4.) One has to estimate the reconstruction error or ac- 
curacy of the reconstructed 3D-structure, which depends on 
the accuracy of the identified structures in 2D, the separation 
angle between the spacecraft, but also on the location, e.g., 
for east-west directed coronal loops the reconstruction error 
is highest close to the loop top. 5.) Eventually we are not 
only interested in the 3D-geometry of loops or plumes, but 
also in physical parameters like density, temperature, plasma 
flow, magnetic field strength etc. Helpful for treating some 
of these problems are coronal magnetic field models extrap- 
olated from photospheric measurements, because observed 
EUV-loops outline the magnetic field. This feature has been 
used for a new method dubbed 'magnetic stereoscopy'. As 
examples we show recent application to active region loops. 



1 Introduction 

The Solar TErrestrial RElations Observatory (STEREO) ob- 
serves for the first time simultaneously the Sun from two 
vantage points (see Kais er et aT] [20 08 ) and allows a three di- 
mensional view of the solar corona. Within this work we aim 
to review stereoscopic methods to reconstruct the 3D corona 
and will concentrate mainly on structures like active region 
loops and polar plumes observed with the STEREO/SECCHI 
instrument package (Sun Earth Connection Coronal and He- 



liospheric Investigation, Howard et al. 2008 1. Before the 
launch of STEREO on October 26, 2006 stereoscopic meth- 
ods have been applied for example to Skylab images by 
IBerton and Sakurail (JT985]l and IBatchelorl (II 994b, and to 



SOHO data by ( |Aschwanden et al.| [19991 12000> . In these 
pre-STEREO cases the authors used the solar rotation and 
took images a few hours to a day apart to reconstruct the 
3D structures under the assumption that all basic features re- 
main stationary within this time. The stationarity assump- 
tion was somewhat relaxed by introducing the concept of 
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dynamic stereoscopy (Asch wanden et aL{|1999| ), which uses 
the a priori information of a coplanar loop shape and a cor- 
responding fitting procedure. In this approach the coronal 
magnetic field -but not necessarily the plasma- is consid- 
ered to be quasi-stationary. The method takes advantage of 
the near-parallelity of adjacent magnetic field lines, even if 
the loop plasma is heated and cools down on much faster 
time scales than the time interval of stereoscopic correlation. 
Such limitations are not necessary anymore after the launch 
of STEREO about two and a half year ago and within this 
paper we would like to give a review on what has been done 
in solar stereoscopy so far, which developments are currently 
under consideration, and planned for the future. 

The key question is how we can derive the 3D geometry 
and physical structure of the solar corona from images ob- 
served with the two STEREO-spacecraft? In (section [2|i we 
describe a step by step guide, which contains the identifi- 
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cation of curvi-linear structures from coronal EUV-images 
in (section [2TT] >, the association of the identified structures 
in both images from different vantage viewpoints in (section 
\2.2) , the geometric 3D stereoscopy (section |2~3| l and estima- 
tion of the 3D-reconstruction error in (section |2.4| i. After 
these steps one has obtained the 3D geometry of, e.g., active 



TRACE 1996_05_19 , 171 A 



region loops or polar plumes and in (section 2.5 1 we outline 



how physical quantities like temperature and density can be 
found. An interesting question, which we address in (section 
|3j, is how well do stereoscopic reconstructed plasma loops 
agree with coronal magnetic field models? Due to the high 
conductivity of the coronal plasma the magnetic field is out- 
lined by the radiating plasma and in principle one has two 
independent data sources about the 3D geometry of coro- 
nal loops, namely stereoscopy and magnetic field extrapo- 
lations from photospheric measurements. In (section [4]) we 
address how in future stereoscopic, tomographic and self- 
consistent modelling approaches could be combined and di- 
minish weaknesses of the individual approaches. 



2 Step by step guide to stereoscopy 

2. 1 Extraction of curvi-linear objects from EUV-images 

A first step for stereoscopy is to extract curve-like structures 
(projections of coronal loops) from observed EUV-images 
such as from the TRACE-image shown in figure[T]a). A prin- 
cipal problem of identifying loops is that the solar corona is 
optically thin and the loops are faint. Visible loops are of- 
ten a superposition of multiple individual loops ( |Schrijver] 
eTaLl[r^9l[2004) l. Figure [T]b) shows the TRACE-image af- 
ter a 7 x 7-boxcar smoothed image was removed from the 
original, which enhances the contrast. |Aschwanden et al. 



(2008a) manually traced 210 loops as shown in panel c). 
While hand-tracing might be suitable to investigate a few 
individual cases, this is not appropriate to study large data 
sets and time sequences. Several automatic feature recog- 
nition methods have been developed. The example shown 
in figure [T] has been used to compare and evaluate five au- 
tomated loop segmentation methods. The output of the five 
codes has been compared with the hand-traced loops from[T] 
c). This comparison revealed large differences. The codes 
identified between 76 and 347 loops. Among the longer and 
more significant loops, the various codes identified between 
19% and 59% of the corresponding 154 hand traced loops 
above this limit. Some codes wrongly identified noise as spu- 
rious short loops. Obviously the state-of-the-art of automatic 
feature recognition techniques is not satisfactory. The indi- 
vidual codes have control parameters, adapted to the image 
signal to noise ratio, resolution and also to the type of objects 
it aims to extract. A prime problem is that time-varying 
background loops and moss prevent an uncontaminated sep- 
aration of loop and background. This makes it difficult to 
trace loop tops and footpoints. The problem becomes more 





Fig. 1. Panel a) shows an original 171 A TRACE image of 19 May 
1998, panel b) a corresponding high-pass-filtered image, where a 
7 x 7-boxcar smoothed image was subtracted from the original. The 
image has been used to compare the output of five feature recogni- 
tion codes with the result of 210 manually traced loop shown in 
panel c). [Original figures are from Aschwa iiden et al.] l |2008a[ > fig- 
ure 1 and 2]. 
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complex for the comparison of EUV-images taken at differ- 
ent wavelengths and temperatures because the background 
is different in each filter. In the current stage the automatic 
feature recognition tools are already useful, but some human 
interaction is usually necessary and the codes can be con- 
sidered as semi-automatic. The automated segmentation of 
polar coronal plumes, which we had during the recent so- 



lar minimum ample opportunity to observe (see Feng et al 



|2009| for details), is more advanced. Even though they are 
fainter in intensity than loops, their shape is more restricted 
so that wavelet and the Hough transform can be employed 



for their detection. Llebaria et al. ( 2002 1 applied the Hough 
transform to a time series (66 hours) of SOHO/LASCO-C2 
observations to generate time intensity diagrams of polar 
plumes. Such automated plume detection can be applied also 
to STEREO-images from two viewpoints (de Patoul et al., in 
preparation). An alternative to stereoscopy is to study polar 



plumes with tomographic methods (see Barbey et al. 2008 



for an application to SOHO/EIT images under the model as- 
sumption of 'polar plumes as stationary objects whose in- 
tensity varies homogeneously with time.'). Two STEREO- 
viewpoints are expected to be most useful for the tomogra- 
phy of plumes when the separation angle is about 60°. 

2.2 Association of objects in both images 

For faint objects like coronal EUV-loops it is not trivial to 
find associated structures in images taken from different van- 
tage points as shown in figure [2] Rodri guez et al.| ( |2009| l de- 
veloped a correlation tracking method which automatically 
matches pixels in both images, which worked well for small 
separation angles between spacecraft but it becomes difficult 
and ambiguous for if the separation angle of the spacecraft 
exceeds about 15°. In particular for large separation angles 
some structures might be visible in one image, but not in 



the other. Aschwanden et al. ( 2008c[ ) applied a forward pro- 
jection with an assumed height range of h = . . . 0.1i? s , 
which was sufficient to find correspondence for 30 traced 
loop, when the separation angle between STEREO-A and B 
was only 7°. The correspondence problem becomes more 
difficult to solve for larger separation angles. 

Other possibilities are the use of a priori assumptions of 
the coronal structures, e.g., fitting to a semi-circular loop 
model ( |Aschwanden et aT [T999), or loop curvature con- 
straints ( |Aschwanden 2005| l" As an alternative one can 
use the fact that the emitting EUV-radiation outlines mag- 
netic field lines due to the high conductivity of the coronal 
plasma. Consequently magnetic field lines should provide a 



reasonable proxy for coronal plasma loops. Wiegelmann and 



|Neukirch| ( |2002| l used the stereoscopic reconstructed loops 
from ( |Aschwanden et al. 1999} and photospheric magne- 
tograms from SOHO/MDI to fit the optimum parameter a 
within the linear force-free field approximation. Projections 
of extrapolated 3D magnetic field lines under different model 
assumptions have been compared with coronal images in 



(|Gary and Alexander! [T999] |Carcedo et al.j |2003| [Regnier 
|and Amain 2004] Wiegelmann et al. 2005 i. The method has 
been extended for a STEREO pair of EUV-images in com- 
bination with linear and nonlinear force-free field models 



Wiegelmann and Inhester (20061 and was dubbed mag 



netic stereoscopy. The idea of magnetic stereoscopy is that 
a number of 3D field line proxies are projected onto the 
EUV-images and compared with the corresponding extracted 
curve-like structures (see section [2~T| . Loops in both images 
which have a minimum distance to the projection of a 3D 
magnetic field line are very likely related to each other. The 
field line proxies where generated from extrapolation mod- 
els but any other method to produce parameterized mean- 
ingful 3D curves would work as well. The extrapola- 
tion models are here just a convenient means to generate 
3D curves the observed loops can be compared with. The 
method has been applied to TRACE-data (taken a day apart) 



by Feng et al. (2007b) and to STEREO/SECCHI by |Fehg 
et al. ( 2007a| l. In both cases the magnetic field has been com 



puted with the linear force free method developed by |See-| 
[haferl ( [19781 from SOHO/MDI. Force-free fields are charac- 
terized by V x B = olB and B-Wa — 0, where B is the mag- 
netic field and a is zero for potential fields and constant in the 
entire space for linear force-free models. For more sophisti- 
cated non-linear force-free field models, a is constant along 
field lines but may vary on different field lines. An advantage 
of using coronal magnetic field models is that they generate 
meaningful and physics-based 3D curves. Note, however, 
that a set of field lines from a linear force-free model con- 
structed with a different a does not constitute a physically 
consistent magnetic field model. Disadvantages are that one 
needs additional observations from ground-based or space- 
borne magnetographs (a third eye, e.g., SOHO or in future 
SDO) and that magnetic modelling and stereoscopy are not 
independent from each other, which is helpful for evaluat- 
ing the consistency of both methods (see also discussion in 
section[3]). 

2.3 Geometric Stereoscopy 

After having solved the feature extraction and association 
problem, the 3D reconstruction is in principle only a geo- 
metric problem as indicated in figure|3]top panel. For stereo- 
scopic reconstruction, it is helpful to introduce a suitable co- 
ordinate system - called epipolar geometry - which reduces 
the original 3D reconstruction problem to a number of 2D 
problems. The positions of the two STEREO-spacecraft and 
any object point define a plane, called epipolar plane, which 
serve as a natural coordinate system for the stereoscopic re- 
construction (see figure [3] bottom panel). The projection 
of the epipolar planes onto the two STEREO image planes 
are called epipolar lines. Per definition all epipolar lines in 
image-A converge at one point (marked as epipole of image 
1 in figure [3] bottom panel) and vice versa. In practice the 
epipolar lines appear almost (but of course not strictly) par- 
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Fig. 2. Top: Contrast enhanced EUVI images from STEREO-B (left) and A (right) of the NOAA AR 0960 on 8 June 2007. The individual 
loops (enumerated white curves) have been extracted by a semi-automatic feature recognition tool as described in Inhe ster et al.| ([2008). 
Please note that equal numbers do not imply a correspondence across the images. Bottom: Some selected coronal loops (yellow) with their 
best fit magnetic field line (red). The separation angle between the spacecraft was 12°. [The top figures have been published originally in 
|Feng et al.| < |2007bt figure 1.] 
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Fig. 3. Top panel: Back projection to reconstruct curve-like ob- 
jects, e.g., coronal loops, from two images. Bottom panel: Epipo- 
lar planes and the corresponding epipolar lines in two STEREO- 
images. [Top and bottom panel of this figure have been originally 
published in Inhester (20061 figure 1 and 2, respectively.] 



allel in images from the STEREO-EUVI (see |Inhesterl [20061 
for quantitative estimations). As the STEREO-spacecraft are 
close to the ecliptic the epipolar lines intersect with the ro- 
tation axis of the Sun (indicated as the z-axis in figure [3] 
bottom panel) and the angle with this axis can be used to 
label the epipolar lines. All points in space, except points 
on the stereo base line (the line which connects both space- 
craft), are lying on a uniquely defined epipolar line. This is 
a helpful constraint as all points, visible on a specific epipo- 
lar line in image-A must lie also on the same epipolar line in 
image-B. In particular, for large separation angles, there is no 
guaranty that a faint coronal structure identified at observer- 
A will be also visible at observer-B. This so called epipo- 
lar constraint (Inhester, 2006) is very helpful because it re- 
duces the 3D-stereoscopic reconstruction to calculating the 
position of each object on it's uniquely defined 2D epipo- 
lar plane. The difference of the positions of an object pro- 
jected along the epipolar lines in the two STEREO-images 
yields the depth information or in case of coronal loops the 
height in the corona. A simple way to tie-point associated 



Fig. 4. Top: How does the uncertainty w of the projected loop 
in both EUV-images affect the error-trapezoid of the reconstructed 
3D-loop? (see text) [Original publishes in (Inhester, 20061 figure 
10], Bottom: Stereoscopic reconstructed 3D loop (yellow) and best 
fit linear force free coronal magnetic field line extrapolated from 
SOHO/MDI magnetograms. [Original published in (Fen g et al.| 
|2007a| l figure 5.] 



loop curves in image A and B is to label along each curve the 
intersections with given epipolar lines and to reconstruct the 
intersection points with the same epipoles line label. This, of 
course, requires that both curves cover more or less the same 
epipolar range, a criterion which can be used to confirm a 
correct association. 

2.4 Estimating the reconstruction error in 3D 

Features in the two STEREO-images can be identified of 
course only within a certain error margin as indicated by w in 
figure|4]top panel. Possible error sources are the finite resolu- 
tion of the instrument as well as uncertainties occurring due 
to extracting features from the EUV-images. The question 
is how do these uncertainties in the 2D-images affect the re- 
construction error of the 3D coronal loop? As a consequence 
of the finite resolution w of the loop projection, the true 3D 
coordinate of a loop point lies nearby in a plane through the 
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reconstructed point which is spanned by the local normals 
of the two projection surfaces, t — (ni x 7i2)/|ni x 7^2! is 
the local curve tangent. Decomposing the uncertainty vector 
into its components n\ and 122 yields an error trapezoid of 
the positional uncertainty has the axes w/ (2cos (a/2)) and 
w/ (2 sin (a/2)), where a is the angle between the projec- 
tion surface normals (see Inhester, 2006] for a mathematical 
derivation and more details). For small a the depth error 
along the mean view direction of the two spacecraft may be 
considerable: 1/ (2sin(a/2)) exceeds 5 for a < 10°. a is 
limited from above by the angle between the two STEREO- 
spacecraft. This upper limit is reached when the loop inter- 
sects an epipolar line normally. This case allows the most 
accurate 3D reconstruction for a given separation angle be- 
tween spacecraft. For features parallel to epipolar lines a is 
zero and the 3D reconstruction error becomes infinite. As 
epipolar lines are almost horizontal in the images this means 
that the highest 3D reconstruction error occurs at the top of 
east-west loops. As a consequence, there are large error bars 
in figure [4]bottom panel, which shown a stereoscopic recon- 
struction of one loop from ( |Feng et aL] |2007a 1 for a separa- 
tion angle of 12° between spacecraft. A small (large) sep- 
aration angle between spacecraft makes it easy (difficult) to 
associate related structures in both images, but for the 3D 
reconstruction error it is the other way around, a small sepa- 
ration angle leads to a large error in 3D. 



An alternative error estimate has been given by Aschwan- 



den et al. ( 2008c 1 for the uncertainty in height. Both error 
estimates show a similar behaviour, however: The smallest 
error is obtained if the loop tangent is normal to the epipo- 
lar planes and the error becomes infinite if the loop tangent 
becomes parallel to the epipolar plane. 

Another problem in finding a unique solution for the 3D 
loop is a possible reconstruction ambiguity. This is a 
problem which can theoretically occur even for a pair of 
correctly identified loop projections, if the footpoints are 
wrongly associated in the image pair. In EUV images, the 
footpoints of loops are sometimes difficult to locate. They 
may be drowned in bundles of other loops or near-surface 
EUV moss. This problem does, however, disappear if one 
can identify the footpoints of loops and requires them to be 
located near the solar surface. It is, however, not always pos- 
sible to identify the loop footpoints and in this case one has 
two possible candidates for the true 3D-solution. The true 
and the false reconstruction intersect at a point, where the 
projected segment is parallel to an epipolar line, e.g., the top 
of east-west coronal loops where the above error estimate 
formally diverges. 

How can ambiguities and errors be limited? One possi- 
bility would be to use additional EUV-images from a third 
viewpoint, e.g., SOHO/EIT or in future SDO/AIA. This pos- 
sibility has (to our knowledge) not been tried out yet. A 
potential problem might be the different resolution of the 
STEREO/SECCHI-EUVI and the SOHO/EUV instruments. 
The error in line-tying is also reduced by making use of the 




Fig. 5. Stereoscopic reconstructed 3D loops and sections of loops 
(yellow) and best fit linear force-free field lines (red) from differ- 
ent viewpoints (top panel: view from STEREO-A, bottom panel: 
Northeast of the Active Region). [Original published in ( |Feng et aT) 
[2007a), figure 3 and 4.] 



fact that the reconstructed loops should resemble field lines 
and hence should be smooth. Smoothing and/or spline fitting 
of the tie point reconstruction should in general reduce the 
reconstruction error. 

Another possibility, which has been already shown to be 
useful for the association problem (section 2.2 1, is to use 



coronal magnetic field models. This possibility, called mag- 
netic stereoscopy was first tested with a model active region 



in ( Wiegelmann and Inhester 2006 ). Different magnetic field 
models, potential, linear and non-linear force-free have been 
used and it was shown that even the use of a poor field model 
(potential fields) was sufficient to resolve the reconstruction 
ambiguity. The extrapolated magnetic field lines provide al- 
ready a proxy for the 3D plasma loop and if ambiguous solu- 
tions occur, the solution closer to this proxy-loop is chosen. 
Figure [4] bottom panel and [5] shows a linear force-free field 
line (extrapolated from SOHO/MDI) in red, which has been 
used for this aim in ( |Feng et aL] |2007a[ ). Extrapolated field 
lines might also serve as a reasonable approximation of the 
plasma loop in regions with a large reconstruction error, e.g., 
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Fig. 6. EUVI-images at different wavelengths from STEREO-A and 
B in the left and right panels, respectively. The dashes white line 
shows a projections of a 3D loop as reconstructed in ( Aschwanden 
|et al.| |2008cl >. [Original published in ( [Aschwanden et al.| [2008b > 
figure 1] 



the loop top or if only parts of the loop are visible in the 
EUV-images. 

2.5 Derive physical quantities 

The geometry of the 3D coronal structures, as seen from dif- 
ferent viewpoints in figure[5]in yellow provide already useful 
information. Feng et al. (2007a} found that most of the recon- 
structed 3D loops cannot be approximated by planar curve 
segments and that most of the loops are not circular. This 
was already known from field modelling. E.g. meaningful 
loop emissions per unit loop length can only be derived from 
EUV images, if the angle between the viewing direction and 



the loop tangent is properly taken into account ( Aschwanden 
|etaLl|2008b) . 

In the following we will concentrate on the computation 
of physical plasma parameters, in particular the electron tem- 
perature and density along the reconstructed 3D loops. Fig- 
ure [6] shows the projection of a 3D loop as reconstructed in 
(| Aschwanden et al.| 12008c]) onto STEREO/SECCHI images 
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Fig. 7. Temperature map of 30 loops. The loops are isothermal 
within the error estimation and the hottest loops are also the small- 
est. [Original published in (Aschwan den et al.| |2008b) > as part of 
figure 9]. 



at different wavelengths taken from STEREO-A and B in the 
left and right panels, respectively. As explained in detail in 
( |Aschwanden et aL~] |2008b[ ) the EUVI-images can be used to 
obtain the electron temperature and density, independent of 
each other for the EUVI-images taken from both STEREO- 
viewpoints, which have been separated by 7° for this study. 
As only a fraction of about 10% of the EUV-radiation is com- 
ing from the loops, one has to remove first the background 
separately in each wavelength, which is a tricky busyness 



and requires model assumptions as explained in detail in ( As 



chwanden et al. 2008b| >. After the well known temperature 
response functions (defined for each wavelength), as calcu- 
lated with the CHIANTI-code, are used to obtain physical 
quantities along the loops. For this aim a local loop-aligned 
coordinate system is introduced and the differential emission 
measure (DEM) is constrained for each loop position with 3 
temperature filters at various positions along the loop with a 
Gaussian function. The DEM is then fitted with the help of 
the EUVI response functions in three wavelength, which pro- 
vides approximations for the temperature, the Gaussian half 
width of DEM and the peak emission measure at the cho- 



sen positions along the loop (see |Aschwanden et aL 2008b 
equations 1-4 for details). Figure [7] shows the temperature 
for 30 plasma loops. The loop-temperatures computed with 
this technique where found isothermal along the loop (within 
an error margin of about 20% caused mainly by background 
substraction uncertainties). The temperature scales with loop 
length with the shortest loops are also the hottest ones. For 
a known loop width (as shown in figure [4j upper panel) the 
three parameters of the DEM-model (temperature, Gaussian 
half width, peak emission measure) can be converted into 
electron densities (see Aschw anden et al . 2008b , section 3.4, 
equations 5-10 for details). With the stereoscopic recon- 
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structed 3D-loop projected onto images from STEREO-A 
and B the temperature and density can be calculated indepen- 
dently from both spacecraft. There have been some system- 
atic differences regarding the results from both spacecraft, 
e.g., the average loop temperature computed from STEREO- 
B was somewhat higher as the one from STEREO-A, but 
the discrepancies have been only a few percent for density 
and temperature estimations. It is still to be investigated 
how consistent the estimations from both spacecraft are for 
larger separation angles. Coronal stereoscopy provides us 
for the first time with the 3D-loop length, temperature and 
density. These quantities can be used to test RTV-scaling 
laws ( |Rosner et al.| [l978) and to go a further step towards a 
self-consistent modelling of the solar corona, which will be 
discussed in section [4] 



3 Stereoscopy and coronal modelling 

As explained in the previous section, coronal magnetic field 
models provide useful information for solving the stereo- 



scopic correspondence problem (see section 2.2 1 and to re- 
move reconstruction ambiguities (section 2.4 1. Until now 
mainly linear force-free models have been used for this aim 
and the method automatically fits the optimum linear force- 
free parameter a for each loop individually to approximate 
its shape as closely as possible, a is here just used as a 
numerical parameter to alter the curve shape. It is there- 
fore not surprising that (Fen g et IT] 2007a| l found a signif- 
icant scatter of a between the loops. We cannot interpret 
the values of a in terms of a linear force-free magnetic field 
model, because the different values are a contradiction to the 
premises of a globally constant a. For a physical meaning- 
ful self-consistent magnetic field model one cannot deter- 
mine the values of a independently for each loop. Unfor- 
tunately, a physical consistent nonlinear model is way more 
demanding, both computationally due to the intrinsic nonlin- 
earity of the model and from an observational point of view 
since nonlinear models require photospheric vector magne- 
tograms as input. Within the last few years a group of scien- 
tists (nonlinear force-free field consortium, chaired by Caro- 
lus Schrijver) has intensively compared and evaluated cor- 
responding computer codes (Schrijver et al. |2006| |2008 



( Metcalf et al.|[2008] l, which showed the codes produce re- 
liable results, when feeded with consistent input data (vec- 
tor magnetograms or a quantity derived from vector mag- 
netograms). In another joint study of the consortium (by 



DeRosa et al. 2009) > the codes have been applied to AR 
10953 and the 3D structure of the magnetic field lines has 
been compared with the 3D geometry of plasma loops (as 



stereoscopically reconstructed in Aschwanden et al. 2008c i. 
Figure [8}, shows the stereo-loops in blue and in red and yel- 
low the magnetic field lines. A major difficulty of this study 
was that the Hinode-SOT/SP vector magnetograms, required 
as input for the magnetic field codes, where available in only 



a very small field of view (dotted area in figure [8]) and has 
only about 10% of the area spanned by the stereo loops. The 
majority of stereo-reconstructed loops were located outside 
of this region, therefore the reconstruction box was signif- 
icantly enlarged beyond the Hinode-area (solid lines in fig- 
ure [8] but only the line of sight component of the photo- 
spheric magnetic field from SOHO/MDI was available in this 
enlarged box. Assumptions regarding the transverse photo- 
spheric field component had to be made in the MDI-area, but 
unfortunately the various magnetic field codes made differ- 
ent use of the assumption in the MDI-area and consequently 
the resulting magnetic field lines differed between the codes. 
The extrapolated magnetic field lines turned out to be also in- 
consistent with the reconstructed STEREO-loops. There was 
an average misalignment angle of 24° between field lines and 
loops. Sand man et al.| ( |2009| ) compared stereoscopic re- 



constructed loops for three active regions observed in April 
and May 2007 with potential field extrapolations and found 
a misalignment angle of about 20° — 40°. In the |DeRosa 



|et al.| ( [2009] l study the nonlinear force-free approach did not 
provide better agreement with the stereo-loops than the much 
simpler to compute potential fields. A main reason seems to 
be the very small field of view, where actually vector mag- 
netogram data have been available and correspondingly the 
codes have not been fed with a consistent input. As visible 
in figure [8] almost no field line/STEREO-loop closes within 
the Hinode field of view. It is therefore necessary to repeat 
such a study with much larger field of view of vector magne- 
togram data (ideally at least the same FOV as spanned by the 
STEREO-loops). 

An additional complication is that the magnetic field vec- 
tor is measured routinely only in the photosphere, where the 
magnetic field is not force-free due to the high plasma j3 
(see Metcalf et al.| 1995| . Consequently photospheric vec- 
tor magnetograms do not provide consistent boundary con- 
ditions for a nonlinear force-free extrapolation. To over- 
come this difficulty, a preprocessing method has been de- 
veloped to remove the non-magnetic forces from the photo- 



spheric vector magnetograms (see Wiegelmann et al. 2006 



|Fuhrmann et al. 2007 for details). These preprocessed 
magnetograms are more chromospheric like and chromo- 
spheric observations, e.g., can be incorporated into the 



preprocessing-algorithm as described in Wiegelmann et al. 



( 2008 ). In principle one could incorporate additional obser- 
vational constraints into the preprocessing routine, for exam- 
ple minimize the angle of stereoscopic reconstructed loops 
(at the footpoints) with the magnetic field vector. This might 
help to better estimate the magnetic field vector in the upper 
chromosphere. Another possibility which might be tried out 
is to add a term which minimizes the angle between the coro- 
nal magnetic field and the reconstructed 3D coronal loops 
in the nonlinear force-free modelling algorithm. One possi- 
bility would be to extend the nonlinear force-free optimiza- 
tion principle Wheatland et al. (2000); Wiegelmann (2004) 
by means of a Lagrangian multiplier £ as suggested in equa- 
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Fig. 8. SOHO/MDI magnetogram of AR 10953 with over-plotted stereoscopic reconstructed loops (see Aschwanden et al. 2008c) in blue 
and extrapolated non-linear force-free coronal magnetic field lines. The solid box depicts the 320 x 256 grid point simulation box in which 
the magnetic field extrapolations have been carried out. The dotted area depicts the field-of-view of Hinode/SOT-SP. Only in this dotted area 
photospheric vector magnetograms have been available. [Original published in (DeRosa et al. , 2009 1 as part of figure 1]. 



tion[T] 

L= f [B- 2 |(V x B) x B| 2 + |V • B| 2 ] d 3 V + 
Jv 

C / (B x S 3D ) 2 d 3 V (1) 
Jv 

The first two terms correspond to the force-free equations 
and the last term measures the angle between the coronal 
magnetic field B and the stereoscopic reconstructed 3D- 
loops S3D, where S3D should contain also an error approxi- 
mation of the stereoscopic reconstruction error This could be 
done by the local length of the vector S3D along the stereo- 
loop, where |S 3 d| = 1 would indicate a small error and 
|S3d| = an infinite error. Locations with |S3d| = 
obviously do not contribute to the functional. The third 
term in (JTJ corresponds to a weighted angle between mag- 
netic field and STEREO-loops. Regions with high magnetic 
field strength and accurate measurement of |S 3 d| contribute 
more to the functional. As such data (an active region ob- 
served simultaneously by both STEREO-spacecraft and a 
large field-of-view vector magnetograph) seem not to be cur- 
rently available (also due to a lack of active regions in recent 
years), corresponding studies have to be postponed until af- 
ter the launch of SDO. The SDO/HMI intmment will pro- 
vide full disc vectormagnetograph, which resolves the lim- 
ited FOV-problem of current instruments. However, the an- 
gle between the spacecraft might then be too large for stere- 
oscopy with STEREO A/B spacecraft. One has to evalu- 
ate whether SDO/AIA and one of the STEREO-spacecraft 
can be used for stereoscopy instead. The result can then 



Where to go in stereoscopy? 
Whe re to g o in corona modeling? 
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Fig. 9. A concept on how stereoscopy could be imbedded in a self- 
consistent coronal modelling approach. 



be compared with nonlinear force-free extrapolations from 
SDO/HMI. This instrument will provide the required larger 
FOV for the field modelling. 



4 Where to go in coronal stereoscopy? 

To summarize the current state of the art of coronal stere- 
oscopy we propose a concept of five steps, namely: feature 
extraction, association, geometric reconstruction, error ap- 
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proximation and physical modelling. Several improvements 
are still possible for some of the steps, e.g., a fully auto- 
matic and reliable feature recognition method, investigations 
on how stereoscopy is still possible with larger separation 
angles between the spacecraft and if we can combine also 
STEREO with other missions, e.g., SDO. Some basic diffi- 
culties, e.g., the large reconstruction error on the top of east- 
west loop is an intrinsic problem of stereoscopy, unless we 
have one or more spacecraft well above or below the ecliptic. 

Despite these principal difficulties coronal stereoscopy 
provides us for the first time some information about the 3D 
geometry and physical quantities of plasma loops. A useful 
concept has been also to combine stereoscopy with magnetic 
modelling of the solar corona, but due to shortage of vec- 
tor magnetogram data these approaches have been mainly 
done with linear force-free methods. With the forthcoming 
full disc vector magnetograph SDO/HMI nonlinear force- 
free magnetic field models are assumed to become avail- 
able on a regular basis. A principal problem of force-free 
magnetic field modelling is that it does not include a self- 
consistent modelling of the coronal plasma. Assume a static 
corona model for example, which balances the Lorentz-force 
with pressure gradient and gravity as j x B = Wp + pVip. 
The term force free means that the Lorentz force vanishes 
and correspondingly the current density j is parallel to the 
magnetic field B. While this approach is well justified for the 
field modelling in the low plasma (3 solar corona, it means 
also that Vp + pVtp = 0, and because the gravity force 
— V?/> is in the radial direction a consequence is that the 
plasma would only be gravitationally stratified, contrary to 
observations. Consequently we cannot describe the solar 
coronal plasma self-consistently within a force-free model. 
Due to the low coronal plasma (3 the force-free assumption is 
well justified to compute the coronal magnetic field structure. 
Variations in the plasma pressure along the field lines are then 
compensated by relatively small Lorentz-forces. For a self- 
consistent modelling of the solar corona we have to take this 
effect into account and use a magneto-hydro-static model in 
lowest order. 

In figure [9] we outline an approach towards a self- 
consistent coronal modelling. First there should be two al- 
ternative routes to derive the 3D geometry of coronal loops, 
either from two EUV-images by stereoscopy or by force-free 
extrapolations from vector magnetograms, e.g., as provided 
in future from SDO/HMI. The resulting 3D field lines and 
3D EUV-loops should be consistent with each other. Unfor- 
tunately this has not been the case for a first comparison (by 
DeRosa et al. 2009 \ as discussed already in section [3] but it 
is assumed (or at least hoped for)that the new instrumentation 
with SDO/HMI and improvements in modelling will lead to 
more consistent results. If consistency between plasma loops 
and field lines has been found, we will have a much more 
reliable magnetic field model than we obtain from extrapo- 
lation alone. One can then derive physical quantities along 



is to use a tomographic approach (as explained in detail in 



Aschwanden et al. 2009 1, which provides the plasma den- 



sity without model assumptions. A sophisticated modelling 
of quantities along the loops can be done (see |Schrijver et all] 



2004 who used scaling laws between the photospheric mag- 



netic field and the heating rate at the footpoints of loops). 
From physical quantities like density and temperature one 
can compute artificial EUV-images and compare them with 
the real STEREO-images. This approach might as well al- 
low to adjust free parameters within the loop modelling ap- 
proach. As pointed out above, modelling of coronal mag- 
netic field and plasma cannot be achieved a selfconsistently 
with a force-free model. Plasma loop modelling will cre- 
ate small gravity and plasma pressure forces, which have to 
be compensated by a Lorentz-forces. Due to the low coro- 
nal plasma j3, the Lorentz force required will be very small 
and the coronal magnetic field structure will deviate only 
marginally from a force-free model. 

Finally a self-consistent model using the magneto-hydro- 
static approach can be computed with the force-free magnetic 
field model and the plasma along the loops as input. Cor- 
responding magneto-hydro-static codes have been developed 
and tested in (Wiegelmann and Neukirch 2006; Wiegelmann 
et al. 2007[> in cartesian and spherical geometry, respectively. 



the loops as outlined in section 2.5 An alternative method 



The resulting consistent static model can be used as input 
to investigate dynamic phenomena for example with time- 
dependent MHD -codes. 

With stereoscopically reconstructed loops, even if they are 
few in number, we will have for the first time an indirect 
though quantitative measurement of the local coronal mag- 
netic field direction. In future magnetic field extrapolation 
codes, this information should be used not only for a com- 
parison but as an additional constraint to lift some of the un- 
certainties which arise from unknown boundary values. For 
a forecast of the evolution and stability of observed active re- 
gions, we probably will need time dependent simulations of 
the corona which incorporate observations by data assimila- 
tion schemes similar to the way they are now used in compu- 
tational meteorology (see Daley 1991| l. An anfolding of the 
solar remote sensing observations to constrain the 3D state 
of the solar corona is vital for these schemes and the forecast 
range will strongly depend on the quality of the 3D recon- 
structions. For this task, the techniques for detecting well 
defined objects in the images and the way they are processed 
into 3D structures needs to be improved. The present state 
of the art of solar stereoscopy can only be a first step in this 
direction. 

A challenging task in stereoscopy is the 3D reconstruction 
of true dynamic phenomena like the initiation of flares and 
CMEs. A key question is which role magnetic reconnection 
plays for such eruptive phenomena. Stereoscopy could help 
us first to derive the quasi-stationary state before an eruption, 
preferably within a sophisticated self-consistent model, e.g., 
magneto-hydro-statics. If we can then also observe the 3D- 
structure of the dynamic phase this would be very helpful for 
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our understanding of such phenomena. One could, in partic- 
ular, use the self-consistent stationary state as input for time 
dependent simulations, e.g., with MHD or Hall-MHD. It is 
notoriously difficult to find transport coefficients, e.g., resis- 
tivity, viscosity, heat tensor etc. for the coronal plasma from 
micro-physics for the coronal plasma, because the required 
kinetic scales are way to small to observe. One could, how- 
ever, try to optimize these transport coefficient with a system- 
atic trial and error approach in order to fit the observations 
best. Such observational based simulations would provide us 
then a rich new world of insights. If the dynamic model simu- 
lations show reasonable agreement with observed quantities, 
we might also have some confidence in other model quan- 
tities, which cannot be observed. This could give insights 
about the physics of magnetic reconnection in the coronal 
plasma. 

A good knowledge of the physics of the solar corona, say 
the 3D structure of loops and plasma along the loops will 
be very helpful, together with high accuracy measurements 
of the photospheric magnetic field vector, to understand the 
interface region between photosphere and corona. A phys- 
ical understanding of this region is important for the coro- 
nal heating problem. Coronal stereoscopy can help here be- 
cause it provides a fair approximation of the corona, but addi- 
tional direct observations of the interface region, as provided 
in near future for example from the small explorer mission 
IRIS are necessary. A modelling approach in this region is 
very challenging, because low and high (3 plasma with sub- 
and supersonic plasma flows exist side by side here. 
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